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Abstract 

We analyze the dynamics of intersecting D3/D3’ brane system overlapping in U-l 
dimensions, in a holographic treatment where N D3-branes are manifested as anti-de- 
Sitter Schwartzschild geometry, and the D3’-brane is treated as a probe. We extract the 
thermodynamic equation of state from the set of embedding solutions, and analyze the 
stability at the perturbative and the non-perturbative level. We review a systematic 
procedure to resolve local instabilities and multi-valuedness in the equations of state 
based on classic ideas of convexity in microcanonical ensumble. We then identify a 
run-away behavior which was not noticed previously for this system. 



1 Introduction 


Recently in [^, a surprising subtelty was identified in a deceptively simple system of inter¬ 
secting D-branes. Consider a system consisting of a D3 and a D3’ brane in type IIB string 
theory, oriented according to 
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and separated by a finite distance in the Xg direction. 

Such a system preserves 8 supercharges. The low energy open string degrees of freedom 
can easily be enumerated as consisting of 

1. An Af = 4 d = 4 U(l) gauge theory living on D3 

2. An Af = 4 d = 4 1/(1)' gauge theory living on D3’ 

3. Two sets of hypermultiplets B and C, arising from Af = 2 d = 4 hypermultiplets, 
dimensionally reduced to d = 1 -|- 1 dimensions and charged as a bi-fundamental under 

u{i) X u{iy 

Explicit coupling between these states was worked out in [^. This system, consisting entirely 
of D3 branes in type IIB string theory, is manifestly self-dual under S'-duality. 

On the hrst pass, there appears to be no obstruction to taking the zero slope limit a' —)■ 0 
as long as one scales the distance separating the D3 branes to be of order 

Axg = a'V (1.2) 

for some V with dimension of mass. This then should give the mass of the B and C helds 
corresponding to the lowest energy 33’ strings. 

As was pointed out in [^, this system exhibits a subtle paradox. A state with a single 
quantum of the B or C helds should exist as a BPS state in the spectrum of the theory, so 
its magnetic dual must also exist as a BPS state in order to be consistent with S'-duality. 
Such a state should arise as a soliton of the held theory at hand. However, the soliton in 
question doesn’t appear to exist; something must therefore be wrong with the assumptions 
being made about the system. 

The resolution proposed by was that the zero slope limit failed to achieve decoupling. 
They then argued that a soliton does exist for a suitably modihed ehective held theory which 
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contains singularities, signaling the need to include additional UV degrees of freedom. By- 
considering full string theory as a UV completion, for instance, the effective held theory can 
be regularized, and the soliton can be constructed as the expected magnetic dual state. 

Following the work of [^, a simple generalization in the brane construction was considered 
in [^. This construction involved also scaling the angle between the D3-branes as 6* ~ a'c, 
introducing a new scale c with dimension of mass squared. This gave rise to a tower of states 
of which B and C are the lightest [^. In that setup, the decoupled theory in the zero slope 
limit is perfectly sensible, and supports the magnetic monopole soliton. One therefore learns 
that one can complete the effective held theory of more economically than by invoking 
full string theory. 

Taking the limit c —>■ oo while keeping V hxed in the construction of essentially 
amounts to recovering the naive zero slope limit considered in . The techniques employed 
in were not particularly ehective for studying this limit, but in [^, we introduced another 
variation in the setup where we replaced the D3 with a stack of N D3’s, so that we have 
as a gauge group U(V) x 17(1)'. This allows us to analyze the system in a strong coupling 
limit where the N D3-branes are replaced by their AdS^ x S 5 dual, and the D3’ is treated 
as a probe. We then considered the magnetic soliton realized as a bion melting into 
the horizon along the lines of (^. It is then possible to see how in the c —)■ 00 limit, the 
magnetically charged bionic soliton delocalizes and decouples as a normalizable state in the 
c —)■ 00 limit. 


The behavior of the D3/D3’ intersection is so counterintuitive that we probably have not 
yet seen the last word regarding this system. One direction which seems potentially interest¬ 
ing is to explore the thermodynamic behavior of this model. An elegant way to approach this 
issue incorprating the effects of interactions is to use gauge gravity correspondince, treating 
the D3’-brane as a probe, along the lines of [^. Working in hnite temperature then amounts 
to studying the embedding of a probe D3’-brane in an AdS-Schwartzschild background. 


Problems of this type where a Dp'-brane probe is embedded into hnite temperature Dp- 
brane geometries have been considered extensively. Most of these works were in the context 
of exploring meson dynamics in holographic QCD [8 15 


It was observed that these brane 


embeddings undergo a phase transition in which they penetrate the black hole horizon. This 
phase transtion was interpreted as “melting” of mesons, which was supported by subsequent 
analysis of the spectrum of huctuations on the probe brane. The behavior near criticality for 
the melting transition also has a rich structure which was noticed even earlier in the context 
of domain walls 16,17. This analysis was carried out in various combinations of Dp-brane 
background and Dp'-brane probes, and the general behavior of the system is dimension 
independent at least in a broad brush perspective. 
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In these class of problems, one generally ennmerates the classical solntions to the eqnation 
of motion corresponding to static brane embedding conhgurations. The embeddings are 
characterized by a control parameter corresponding to qnark mass, and have a dehnite 
order parameter corresponding to the quark condensate. The static solutions constrain the 
equation of state relating the quark condensate to quark mass. Treating these parameters 
as thermodynamic quantities, one can explore issues such as thermodynamic stability and 
hydrodynamic limits. Indeed, these system generally exhibits instabilities and multi-valued 


equation of states as were observed, for instance, in 18 


The case of D3’-brane probe in the background of D3 branes oriented according to (1.1), 


however, is somewhat special. This is the case that was singled out in Reference [21] of 10 


There are two concrete senses in which the D3/D3’ system stands out. The brane embedding 
is characterized by a scalar on the D3’-brane world volume which happens to saturate the 
Breitenlohner-Freedman > — 1. (^. As a result, the operators corresponding to these 
helds have scaling dimension 


^ d =F V(P + 4777-2 
A = -= 1 


(1.3) 


which is degenerate for d = 2 and = —1. This gives rise to logarithmic factors in the 
scaling of the embedding solution near the boundary, and exchange the roles of control and 
order parameters in the holographic dictionary as we will ellaborate further below. Also, the 
holographic renormalization procedure requires introducing an anomalous scale which can 
affect certain physical observables (22). 


The thermodynamics of D3/D3’ system have been analyzed previously 23 although 
in that threatment, the control parameter of the embedding was treated as being hxed 
in dehning the ensumble. The analysis of 23 was primarily focused on establishing the 
existence of a robust zero mode in the longitudional electric fluctuations and its implication 
for charge transport behaviors. 

In this article, we re-examine the thermodynamics of D3/D3’ system with emphasis on 
understanding the thermodynamic stability issue of the quark condensate order parameter. 
We will explore the stability both at the perturbative and the non-perturbative level, and ar¬ 
gue that the phase diagram of the system looks somewhat different than what was suggested 
For now, we will focus primarily on zero charge embeddings. The extention of this 


m 


23 


analysis including charges and chemical potential will be reported in a separate publication. 
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Figure 1: A schematic illustration of a flat brane embedding. This embedding is static in 
the absence of the black hole. We are interested in how this embedding is deformed when 
the black hole is introduced. 

2 D3’-brane probing finite temperature anti de Sitter space 


In this section, we will review the analysis of embedding D3’-branes in the near horizon 
geometry of N D3-branes at hnite temperature. Similar analysis can be found extensively 
in the literature, but it is useful to formulate it here to make the notation and conventions 
explicit. 

We begin by writing the supergravity solution corresponding to a stack of N D3-branes 
at hnite temperature in type IIB supergravity 


ds^ = + dx^) + + r^dnl) 


with 


and 


d4 

// = ! + —, = ATig,Na'^ = \a'^ , 




By standard arguments, the temperature T and the horizon radius Vg are related by 

T= ^ 


TlTg ttaA 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 


We will take the decoupling limit by sending a' —)• 0 keeping U = r/a' and T hxed. This 
will also keep Ug = Xg/a' hxed and hnite. 

At zero temperature, brane conhguration oriented as (1.1) and illustrated in hgure[^is a 
consistent static solution. We are interested in how such an embedding is deformed when T 
is no longer zero. 
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We will parameterize the six dimensions transverse to the N D3-branes as 


X4 

= r sin 01 sin 02 sin 63 sin 64 cos 05 

(2.5) 

X5 

= r sin 01 sin 02 sin 63 sin 04 sin 05 

(2.6) 

xe 

= r sin 01 sin 62 sin 63 cos 04 

(2.7) 

X7 

= r sin 01 sin 02 cos 63 

(2.8) 

Xs 

= r sin 01 cos 02 

(2.9) 

Xq 

= r cos 01 . 

(2.10) 


We can then treat 

t, z = X3, r, 0 = 6^5 (2.11) 

as the world volume coordinates in static gauge for the embedding, and denote 

d^a = dt dz dr d(j) r (2-12) 


and treat 6 i, 62 , Os, and 64 as parameterizing the embedding. We will restrict our attention 
to conhgurations where the gauge held on the world volume of D3’ is trivial at hrst. The 
system is invariant under S'0(4) rotational invariance acting on (xg, xy, Xg, Xg), and it turns 
out to be dynamically consistent to set all but one of the four coordinates to zero. This is 
equivalent to setting 

02 = 03 = 04 = ^ (2.13) 

and treating 0 = 7r/2 — 0i as the only relevant held variable. Restricting to embeddings 
which are invariant under translation in t and z directions, the resulting ehective action is 

^ = {o /2 [ cos(0(r))^l + r2/(r)0'(r)2 . (2.14) 

[ZTT) Ol Qs J 

Note in particular that the dependence on warp factor H dropped out completely from the 
act ion 0 

In the zero slope limit, one sees that a' scales out of the action 


Sdbi — 


1 


d^E cos(0(f/)) 




9'{Uy 


(2.15) 


with 


= dt dz dUd^U . 


(2.16) 


^The warp factor is nonetheless relevant for the formula relating temperature to horizon radius (2.4). The 
warp factor also enters in the computation of quasinormal modes in appendix [Bj 
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For the purpose of finding the embeddings which extremizes this action, it is convenient 
to scale Us out by dehning 

(2.17) 


u = 


U 

TTs 


so that the action takes the form 

t/2 '• 




where L is the volume of the Xi coordinate and T. 


d X du u cos(6'(m)) 4 /1 + -r 6 '{uY 




(2.18) 


The task at hand now is to solve the equation of motion obtained by varying (2.18) which 


is a non-linear second order differential equation for 9{u) for u taking values in the range 
1 < M < oo. For large values of u, 9{u) approaches zero and one can show that the solution 
can be parameterized in the form 


9{u) = 


clog('u) 


u 


m 

H- 

u 


(2.19) 


where we denote dimensionless integration constant^ c and m following the convention 
of ( 23 ]. Near u = 1, we impose the regularity condition which constrains c as a function of 


m. 


The actual solutions satisfying these boundary conditions have to be obtained numeri¬ 
cally. The general feature of the solutions we obtained is illustrated in hgure There are 
two different class of solutions depending on whether the brane penetrates the black hole 
horizon or not. The ones which do not, known as Minkowski embeddings, are illustrated 
in red. The ones which do, known as black hole embeddings, are illustrated in blue. Each 
of these solutions have a definite value of c and m. The set of (c, m) computed for these 
solutions are illustrated in figure 


There are a number of features that are worth noting in hgures and First, m is 
not single valued as a function of c. Closely related is the fact that there is a maximum 
value of c where dc/dm = 0. Also, note that the embedding exhibits a self similar critical 
structure when the Mikowski and black hole embeddings meet. This was a feature originally 
noted 16,17 and further ellaborated in the context of Dp/Dg system in 10,11 • It signals 


that there is a hrst order “meson melting” phase transition near the self-similar critical point. 
In this article, we will have more to say about the critical behavior at dc/dm = 0 than at 
the self similar point. 


^This c is unrelated to the c parameterizing the tilt of D3 relative to D3’ in the notation of [^. 
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Figure 2: Black hole embedding illustrated in blue and Minkowski embeddedings illustrated 
in red. 
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Figure 3: {c,m) for Minkowski (red) and black hole (blue) embeddings of D3’-probe. Each 
dot corresponds to numerical solutions we found with initial conditions for 0{U) specified 
either at the horizon U = Ug for the black hole embeddings, or by fixing 6 '{U) at 0{U) = 7r/2 
for the Minkowski embeddings. 
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3 Thermodynamics and Holography of D3/D3’ system 


In this section, we elaborate on the thermodynamic and holographic interpretations of em¬ 
bedding solutions illustrated in hgures and 

First, it should be noted that all embeddings illustrated in hguresandare asymptoti¬ 
cally AdSs X Si- The world volume degrees of freedom on D3’-brane will have a holographic 
interpretation as a held theory in 1-1-1 dimensions. 

The embedding held 0{U) will be associated to an operator of dimension A = 1 to be 
associated with quark bilinar via the standard holographic dicationary which needs to 
be stated with some care because A is degenerate. Let us ellaborate on this matter further. 

To interpret the system holographically, it is awkard to scale out Us since the holograh- 
pic dictionary should be independent of Ug. Let us therefore parameterize the asymptotic 
behavior of the 0{U) in the form 


9{U) = 


Clog{U/U, 

u 


M 

+ T7 


(3.1) 


where once again adapting the notation of 23 , C and M have dimension of mass, whereas 
t/* is an arbitrary scale with dimension of mass which one must introduce in order to make 
sense of the argument of the logarithm. An astute reader should notice at this stage that 
there is some ambiguity in how M is dehned since changing 17* has the effect of shifting 
M. It would therefore be essential to understand if and how t/* affects physical observables 
(or not.) Several related issues will arise in the discussion below and we will be doing due 
diligence to track these issues. 

We begin by recalling the standard formulation of the holograhpic dictionary (see e.g. 
for a review) that 

Zbuik[C{x)] = (3.2) 


where the left hand side of the equality describes a path integral for bulk helds such as 6 {U) 
carried out in such a way that 9{U) asymptotes to 


9{U) ~ 


C{x) log(17/17*) 
U 


(3.3) 


as the boundary is approached by taking f/ —)■ oo. This path integral and boundary condition 
as formulated is independent of f/*. One might formulate the right hand side for the 9{U) 
held to take the form 




(3.4) 





where Sdbi is the action given in (2.15), and the boundary condition for 6 {U) is referenced 
implicitly in the specihcation of the measure. One can apply the saddle point approximation 
to identify the domimant contribution to this path integral, which simply amounts to evalu¬ 
ating the action for the solution to the equations of motion enumerated in £gures[^and|^ As 
is typical in these computation, however, the action formally diverges, and a renormalization 
is required to dehne the bulk side of the correspondence unambiguously. For the D3/D3’ 
system, this was worked out explicitly in (6.2)-(6.4) of 22 . We are simply instructed to add 
the holographic renormalization counter-term 


ScT — 


With this counter-term included. 


d^xU^ - 1- 


\og{U/UcT] 


9{uy 


(3.5) 


U=Uuv 


Zbuik[C{x)] = j [ 110 ( 17 , ( 3 . 6 ) 

where Ujjv is the ultra-violet cut-off scale, whereas Uct is a new scale that is required in 
order to make the logarithm appearing in the counter-term make sense. This expression 
is hnite in the Uuv oo limit. However, the dependence on Uct which characterizes the 
renormalization scheme survives and should a priori be treated as independent of f/* and 
Ug. A useful feature to isolate in the holographic dictionary is the prescription to extract 
the expectation value of the operator dual to 9{U). This is to be derived by varying the 
logarithm of Z[C{x)\ with respect to C{x). By manipulating Sdbi + 5'ct, one hnds that 

{0(i)> = - C(x) \og(U.IUcT)) ■ (3.7) 

{2'Kygs 


It is also useful to infer the values of (C, M) for the solutions enumerated in figures 
andj^which are parameterized in terms of (c,m). By simply relating (2.19) to (3.1), we hnd 
that 

C = UsC, M = Us{m-c\og{Us/Uy) . (3.8) 


In terms of (c, m), we can write 

{0{x)) = -clog{Ug/UcT)) (3.9) 

whose signihcance is the fact that the dependence on f/* has dropped out. However, the 
dependence on Ug and Uct remains. 

At this point, we can also compute the free energy by evaluating the action with time 
compactihed on a circle of radius 1/2tiT 

G{C,T) = -T\og{Z[C]) (3.10) 
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Figure 4: M{C,T) and G{C,T) as a function of C for \og{Us/UcT) = —2 (top) and 
\og{Us/UcT) = 2 (bottom) 


or by computing 

G(C,T) = [ dC'(M(C)-C'\og(U./UcT)) 

Gt^YQs Jo 

LiP P 

= ~ (0 Y / c?c'(m(c) - c'log(t/,/[/cr)) (3.11) 

GT^rgs Jo 

which gives an equivalent independent result. Since is essentially unphysical, it is 
conveninet to set = Uct so that 

for the remainder of this paper. 


We are also now in the position to display thermodynamic data such as the equation of 
state M{C, T) and the free energy G(C, T) for various fixed values of T/Uct- Few examples 
are illustrated in figures 


It is also straight forward to infer quantity such as 


S{G,T) 


dG 

W 


LUs 

{27rYgs 


m{c)dc! 



(3.13) 
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Figure 5: S{C,T) as a function of T for fixed C. This curve turns out to not depend on 
UcT- 



which happens to be independent of Uqt and plot it as a function of T with C hxed. We 
illustrate that in figure This is essentially equivalent to what is illustrated in hgure 3.d 
of 23 except that we do not hnd any evidence of the dotted part of their graph. This view 
is supported also from the structure of equation of state illustrated in hgure which has 
exactly two, not three, branches as c approaches zero. We will ellaborate further on this 
point below. 


4 Physical interpretation of D3/D3’ thermodynamics 


In the previous section, we outlined the general features of the D3/D3’ system which can be 
presented in the thermodynamic context. The control and order parameters, {C,M) play a 
role very similar to mechanical parameters (P, V), magnetic parameters {H, M), etc., up to 
a sign which arises from conventions which are set for historical reasons. 


There are two notable features about the equation of state and the subsequent thermo¬ 
dynamics summarized in the previous section. One is the fact that the scale Uct affects 
the equation of state which is physically observable. The other is the fact that the equation 
of state exhibits multi-valued-ness and regions of instability. This latter issue is somewhat 
familiar from previous consideration of black hole thermodynamics 18 21 . It is generally 
stated that the van der Waals model of liquid-gas phase transition is a prototype for under¬ 
standing these issue. Nonetheless, in the context of van der Waals model, it was the order 
parameter as a function of the control parameter, V (P) that was multi-valued, whereas in 
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Figure 6: F{M) for log(f/s/f/cr) taking values —2 (left), 0 (center), and 2 (right). For all of 
these cases, F{M) asymptotes to a constant value 0 as M is taken to infinity. 


the case of D3/D3’ system, it is exactly the opposite]^ Another apparent paradox stems 
from the fact that the sucsceptiblity 


_ dC 
^^~dM 


(4.1) 


characterizing this system is explicitly dependent on Uct whereas dymamical features such 
as the poles of quasi-normal modes at fixed control parameter C are manifestly independent 
of Uct- The goal of this section is to clarify these issues. 


Let us begin by recalling the classical thermodynamic perspective on stability. A useful 
quantity to consider is the effective potential of the order parameter obtained by Legendre 
transforming the free energy 


F(M) = G{C) 




-CM 


G'(C)=-M 


One can also compute F{M) in terms of the equation of state 


fM 


F(M) = 




dM'C(M') 


( 4 . 2 ) 


( 4 . 3 ) 


Plotted as a function of M with hxed values of \og{Us/UcT), they take the form illustrated 
in figure]^ For all of these figures, M —)■ oo corresponds to the Minkowski branch, and F{M) 
approaches a constant, reflecting the fact that the area under the curve C{M) is finite. 


The free energy F{M) is an important physical quantity characterizing the effective 
thermodynamic behavior of the order parameter M. If a system with free energy F{M) is 
brought to contact with a reservoir with which the M-ness is freely exchanged, the system 
will achieve equilibrium when M minimizes the potential 


T,f(M) = F(M) - 


^For an illuminating discussion of this very issue, see 


26 


( 4 . 4 ) 
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where Cext is the control parameter conjugate to M of the reservoir. One of course recognizes 


G{Cext) = F{M) - 


CpxfM 


(4.5) 


F'(M)-Cp:,t=0 


{2TiYgs 

as the conjugate Gibbs Free energy when the M that minimizes Feff{M) is substituted into 
Feff(M). From that point of view, it is natural to associate the susceptibility 


(27r)^5fs , dC 

y = f— L^F''{M) = - 

L ^ ’ dM 


(4.6) 


as parameterizing the stability of the system. The susceptibility xe is dependent on Uct- 

Let us now take a closer look at the form of F{M) illustrated in figure]^ and make several 
observations. 


1. F{M) is multi-valued over some range of M. 

2 . The susceptibility at M = 0 (as well as other values of M) changes as Us /Uct is 
varied, and can get negative, signalling an instability, for instance for large positve 
values of \og{Us/UcT)- 

3. The effective action Feff{M) is not concave everywhere and does not have global stable 
minima when is non-vanishing. 

Let us address each of these observations more carefully. 


4.1 Multi-valuedness of F{M) 

This issue is not too serious. The fact that there are multiple branches for some hxed 
value of M (and T) simply reflects the fact that there are multiple thermodynamic states 
corresponding to these order and hxed parameters. However, in thermodynamics, one focuses 
on the dominant state in the ensumble, which is the one with the lowest free energy. So in 
reading hgurej^ one should simply trace the branch with smallest F{M) for any hxed value 
of M, regardless of the discontinuities that might result. 


4.2 Susceptibility and its dependence on Uct 

This is an extremely important yet subtle issue. Taken at face value, it implies that the 
susceptibility and therefore the thermodynamic stability depends on Uct- For example, in 
hgure 1^ we see for \og{Us/UcT) = 2 that F{M) is concave down indicating instability at 
C = M = 0. 
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On the other hand, it is generally established that the thermodynamic stability can be 
inferred from the presence or absence of poles of quasi-normal modes in the upper half of 
the complex uj plane 27 . The quasi-normal mode analysis, however, does not depend on 


holographic renormalization counter-term. As such, it would appear that thermodynamic 
stability is independent of Vct- But this is in direct contradiction with what we stated in 
the previous paragraph. 


The eventual resolution of this apparent tension can be understood as susceptibility 
being dependent on Uqt but not the stability. But there are number of subtleties involved 
in arriving at this conclusion which we will describe in this subsection. 


The issue boils down to mapping out the allowed range of values in enumerating the 
renormalization schemes which XJqt is parameterizing. On the other hand, Xe{M = 0 ) 
depends on Uqt and as such can also be considered as parameterizing the renormalizaiton 
schemes. The susceptibiltiy Xe{M = 0) however is a quantity that is easy to measure. Uct, 
on the other hand, is an abstract quantity appearing in the counter-term which can only be 
inferred by measuring some physical quantity (such as Xo{M = 0)) and using its relationship 
to Uct- The choice to parameterize the renormalization scheme with at M = 0 is an 
arbitrary choice. Any other M can be used as a reference. The situation is analogous to the 
relation between renormalized coupling in MS scheme and physical coupling inferred from 
scattering at some definite energy. The former is the analogue of Uct whereas the latter is 
the analogue of Xe{M = 0). 

The question then is what constitutes the appropriate range of parameters to enumerate 
distinct renormalization scheme. Should it be 0 < Uct S oo, or —oo < Xe{M = 0) < oo7 
To the extent that Xe{M = 0) is the physical parameter, it would seem natural to treat 
the latter as parameterizing physically distinct renormalization schemes. We will adopt that 
point of view in this paper. There is a possibliy that extrapolation beyond infinite Xe{M = 0) 
would admit interpretation along the lines of dualities where one extrapolates beyond inhnite 
coupling, but we will not pursue that possiblity in this paper. 

This implies however that the behavior illustrated in figure and for \og{Us/UcT) = 2 
where the Xe{.M = 0) is negative signalling instability corresponds to pushing xe{.M = 0) 
beyond inhnity and should be excluded from our analysis. At M = 0, one can explore the 
full range of 0 < Xe{M = 0) < cx) by letting Us/U ct vary. Susceptibility at M = 0 for this 
system is always positive. 


This however does not imply that the system is always stable or that the stability analysis 
is completely unrelated to the quasi-normal mode analysis. To see this, suppose we set the 
temperature Ug = Uct so that the equation of state is given by what is illustrated in hgure 
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L(U,f 



Figure 7: = F{M) - Ce^tM with C/Us = 0.1. 


For every 0 < C < Cmax, there is a subleading branch of solutions where 


dC 

dM 


< 0 . 


(4.7) 


One can see the same thing by looking at unstable stationary point in the potential illustrated 
in hgure which is appropriately tilted by the inclusion of —CextM term as is illustrated 
in fugre One exepcts to hnd a corresponding pathology in the spectrum of quasi-normal 
modes for the fluctations around such unstable background solutions, but how can one 
understand its appearance short of doing an explicit computation? 


One way to see that a cross-over into unstable behavior is taking place is to focus on the 
stat^at C = Cmax where xe = dC/dM = 0. Recall that for every point on hgure there is 
a corresponding solution 6 {U) which extremizes the action (2.15). If we parameterize these 
solutions for hxed value of M by 9m{U), then at M = Merit corresponding to O = Cmax, h 
follows that 


iiU) = 


(4.8) 


M=Mcrit 


‘ip{U) is a gap-less quasi-normal mode with u = k = 0. The fact that a gapless mode is 
appearing precisely when the susceptibility Xe = ^ strongly suggests that an unstable pole 
would appear upon continuing to the branch where the susceptibility is negative. This cross¬ 
over across = 0 is different from crossin over xe = oo which we discussed earlier in the 
context of placing a bound on Uqt- 


"^Note that the position Cmax is independent of Uct since a change of Uct can only shift Xg ^ by a finite 
amount. 
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By a similar token, for sufficiently large negative value of log(t/s/t/cT), we encounter a 
point in the unstable branch where xe reaches negative inhnity. This can be seen in hgure 
1^ where the tangent of M{C) is horizontal. For \og{Us/UcT) = ~2, this point corresponds 
to the switch back point in the black hole embedding branch illustrated on the left most 
hgure otg This issue is somewhat academic, however, since the switch back point is already 
subdominant in the saddle point approximation as can be seen in hgure 


The critical point (C, M) = {Cmax, Merit), on the other hand, is a dominant saddle point 
and as such the resulting critical behavior giving rise to a new channel for dissipating energy 
in the hydrodynamics limit near that point is a real physical feature of this model, at least 
at the perturbative level. This critical point can also be seen to correspond to the point in 
hgure 1^ where T = takes on a minimal value. The entropy S plotted as a function of T 
in hgure is double valued. From the thermodynamic point of view, the dominant branch, 
however, is naturally the one with greater entropy. So, the black hole embedding dominates, 
and the Minkowski embedding is the subdominant one. For T > T^m, the system naively 
seems to be perfectly stable and well behaved. At T = Tmin, there is a critical behavior. 
The phase strucure implied by these features are consistent with the d = 0 slice of the phase 
digram illustrated in hgure 4 of 23 . This however raises one obvious question. If at T = Tmin 
we encounter a crticial behavior, where does the system equilibriate to for T < Tmin- In 
order to address this issue, we need to go beyond the scope of perturbative stability analysis, 
and consider the global issues. We will discuss that issue in the next subsection. 


The apparent quantitative mismatch in the dependence of Uct between susceptibility 
and quasi-normal mode specturm can also be seen in the computation of correlators in the 
real-time formalism. The retarded Greens function is computed using the prescription given 
in equation (3.15) 

G\k) = -2:F{k,z)l, = V^g^^U{z)dMz) (4.9) 

for a suitably normalized ingong wave f{z), where z ~ 1/u. However, strictly speaking, T 
is divergent for our model and a counter-term is needed to render this expression hnite. 


One can understand the origin of this divergence as arising from computing the two ponit 
function of operator whose dimension is A so that at short distance, it scales as 

G{k) = {O(x)O{0)) ~ ^ (4.10) 

which then in momentum space takes the form 

Ol.k) ~ ^ (4.11) 

for large k. The issue arises when d — 2A < 0 so that this correlation funciton do not decay 
at large k. This is the case in our example because d — 2A = 0. So strictly speaking, one 
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expects the two point function to scale for large k as 

G{k) ~ log(fcVh^) (4.12) 

for some scale /r. We are however interested in the small k behavior when the system is at 
hnite temperature. 

The two point function should then admit a spectral decomposition 

G'(fc) = / (4.13) 

where by power counting, we know that p{s) must asymptote to a constant at large s. The 
integral over s, however, does not converge and must be regulated, for instance, by adding 
a term 

G(k) = J ds =ao + ^2^^ + a4k^ + ■■■ (4.14) 

The term added is a contact term in that it is independent of A;. It only affects the Oq 
term in the small momentum expansion of G{k). The scale Uct and p arises from these 
considerations, which does affect the two point function and therefore the susceptibility, but 
does not affect the pole structure of G(A;)|^ Nonetheless, both the susceptibility and the 
quasi-normal mode spectrum knows when the system is perturbatively unstable, and exhibits 
the appropriate symptoms. 

4.3 Non-perturbative stability of the D3/D3’ system 

In this section, we will discuss the subject of how the unstable states relaxes to the true 
equilibrium state. This issue requires consideration beyond the perturbative analysis, but 
the subject is not an unfamiliar one. The same issue arises in the phase structure of liquid- 
gas transitions in the van der Waals model. Let us see how that applies to the D3/D3’ 
system under consideration. 

It is a fundamental fact of statistical mechanics that the set of accessible states param¬ 
eterized by the order parameters form a convex set. A nice historical review of this basic 
notion can be found in [^. If one is working at hxed temperature T, in a system with 
a single order parameter M, the region bounded by the curve F{M) must be convex, or 
equivalently, F"{M) > 0. But sometimes, as is the case here, by working out the equation of 
states for what one believes is the dominant thermodynamic conhguration, one hnds regions 
where F"{M) is negative, signalling an instability. For our system, this can be seen very 
explicitly in hgurej^ where F{M) is concave down for large values of M. Equivalently, one 
can see regions where dG/dM is negative in the equation state illustrated in £gure[^and|^ 

^The issue of subtle contribution from contact terms was also discussed in (2.12)-(2.14) of [^. 
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Figure 8: Typical equation of state P{V) for various fixed T illustrating thermodynamic 
instability, coexistence of states, and the critical point. As T is varied, regions of instability 
where — dPldV\j, < 0 appears and is shaded in red. This causes the system to undergo a 
hrst order phase transition where the equation of motion is modihed by the ruled coexistence 
line, illustrated in black. The coexistence region is generally larger than the unstable region. 
So the local instability in the red shaded region affects the equilibrium conhguration of all 
points in the (P, V) region inside the blue shaded region. 

The standard procedure when this happens is to observe that the system can acheve 
a state with lower free energy for fixed M by being in a hetrogeneous co-existence state. 
A co-existence of two state will give rise to a configuration with (M, P(M)) interpolating 
between the pair of states. As a result, the maximal extension of the space of states al¬ 
lowed by considering co-existence states is precisely the convexification of F{M) achieved by 
supplementing F{M) with “ruled surfaces” in the terminology of [^ . 

This is exactly what happens to the equation of state in van der Waals system. We 
illustrate the standard diagram displaying a collection of isothermal P{V) curve in figure 

The point to note is the fact that 1) there are regions where —dP/dV < 0 signalling 
instability, and 2) that this gives to a modihed equation of state by allowing coexistence 
states. It should also be emphasized that 3) the region of phase diagram modihed by allowing 
coexistence regions (shaded blue region in hgurej^ are strictly greater than the region where 
the system exhibts apparent instability (shaded red region in hgurej^. 

How do these ideas apply to the D3/D3’ system? From the equation of state illustrated 
in hgures and 1^ where C{M) is rapidly approaching zero as M is increased, it follows that 
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F{M) must approach a constant value (set to zero for convenience in figure]^) As soon as 
this potential is slighlty tilted in response to non-vanishing external control parameter C^xt 
as is illustrated in hgure the displaced local minima is no longer a global minima and the 
system is susceptible to decay via tunneling into a run away behavior towards large M. If 
Cext is taken to be larger than Cmax, even the local minima disappears and the potential 
does not have any stationary points. (The situation is a little different when the some net 
charge is introduced to the D3’-brane world volume. We will ellaborate further on this case 
in a separate publication.) 

From the point of view of convexifying F{M), we see that the result is to say 

F{M) = F{M = 0) = constant (4-15) 


which is even more susceptible to runaway behavior when C is non-vanishing. 


What appears to be happening is that the d = 0 slice of the phase diagram illustrated 
in hgure 4 of 23 degenerated completely and that at hnite C, the system is completely 
unstable at the non-perturbative level. 


What this means, presumably, is that the classical treatment is predicting its own ther¬ 
modynamic demise and that some conhguration not presently accounted for will modify the 
equation of state to provide the stable, equilibirium state for this system. Such corrections, 
however, must arise from string or quantum corrections and is expected to scale non-trivially 
with respect to gs and N. It is possible that such correction would also give rise to the branch 
drawn with dotted lines in hgure 3.d of (2^ . 

Alternatively, the D3/D3’ system is intrinsically unstable. At this moment, we do not 
have any reason to rule out that possiblity. 


In this article, we took the boundary of AdS^ x to be hat and inhnite in volume. It 
would be interesting to repeat this analysis treating the boundary to be a S '3 of hnite size. 
In that case, the entire system undergoes a Hawking-Page transition 30 , giving rise to a 
qualitatively new behavior also for the D3’-brane embedding. It is somewhat unexpected, 
however, for thermodynamic stability of a system to rely on hnite volume issue. In any 
case, it would be very interesting to understand all the diherent ways in which the run-away 
behavior seen in this system can be stablized. 


5 Discussions 

In this article, we analyzed the embeddings of D3’-brane probe in Schwarzschild AdS^ ge¬ 
ometry and studied their thermodynamic interpretations, with emphasis on order parameter 
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and control parameter dual to the brane probe embedding psuedo scalar field 9{U,x). We 
described the subtle relationship between the anomalous scale Uct introduced in the holo¬ 
graphic renormalizatoin procedure, thermal susceptibitliy xe, and the spectrum of quasi¬ 
normal modes. 


The resulting analysis of the theromdynamic stability revealed that while the system is 
perturbatively stable for some range of parameters, it is unstable at the non-perturbative 
level almost everywhere. The full implication of this instability is not completely clear to us 
at the moment. It should be noted that our consideration was limited to treating the branes 
in the probe approximation. Perhaps gravitational back reaction will stabilize the system, 
although properly addressing this issue is a tall order for an intersecting brane system. It 
is also possible that a satisfactory resolution will require going beyond the semi-classical 
treatment of these systems. 


The original goal of this study was to find some relation between subtle features of the 
thermodynamics of D3/D3’ brane system to the subtle features discussed recently by Mintun, 
Polchinski, and Sun [^. At the moment, the only connection we see is the fact that some of 
the pathologies are due to low co-dimension physics in both instances. We hope to provide 
more insight into this issue in the future. 


Finally, let us comment that the intersecting D3/D3’ system could also prove to be 
useful as a probe of the region behind the horizon in the context of entanglement entropy 
where the degrees of freedom crossing the horizon is in the open string sector IsTl. Just like 


in other attempts to probe behind the horizon such as 32,33 , we expect the sensitivity of 
probe dyanmics to behind the horizon physics to be highly suppressed. Nonetheless, perhaps 
something can be gained by using an open string probe instead of a closed string probe. 
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A More on embeddings 


In this appendix, we offer an alternative parameterization of the equation of motion (2.15) 
where we map the problem to dynamics of a particle rolling down a potential. This formula- 


20 






tion is useful for visualizing the solutions and for providing assurance that all solutions are 
accounted. 


The procedure to convert (2.15) into the classical potential problem form is essentially 


the same steps one takes to convert the Nambu-Goto action in to the Polyakov action form. 
This can be done by introducing an auxiliary world volume parameter r and a Lagrange 


multiplier A(r) and re-writing the effectively one dimensional form of (2.15) as 


dr-G(r)cos(0(r)) [A(r) + {U{Tf - U^U{t) 2)d'(r)2) + A(r)] . (A.l) 


Solving for the Lagrange multiplier and setting t = U recovers (2.15). On the other hand, 
imposing as the gauge condition 


[/2 _ t /4^-2 

A(r) = •- ^ - Ucose , 


U! 


setting 


= f/f cosh s , 


and rescaling r = UgCr will scale out Ug, leads to 


L = [ da -I- -I- - sinh^ s cos^ 9 

\ 2 2 2 


(A. 2 ) 


(A.3) 


(A.4) 


where s, d, and a are dimensionless, and dot denotes derivative with respect to a. In the 


form (A.4) the problem is essentially that of a particle whose positions are parameterized by 


1 < s < oo and — 7 r /2 < 6 < n/2 rolling along a potential 


V(s, 9) = — - sinh^ s cos^ 9 . 


(A.5) 


Reparameterization invariance further implies that the solution should correspond to trajec¬ 
tory with vanishing Hamiltonian 

n = 0 . (A. 6 ) 

A large class of solutions arises as a trajectory of a particle rolling up, turning around, 
and rolling back down the potential in (s, 9) coordinates (A.5) which we illustrate in hgurej^ 


These are the worm hole embeddings in the terminology of 16 OT . In the original (X 4 5, Xg) 


cooridinates, they look like an embedding illustrated in hgure 10 In the large X region. 


these correspond to having both a brane and an anti-brane along the lines of [6,34 and 
not the solutions we are looking for. 


are 


The only other possibility is for the trajectory in the (s, 9) plane to hit the boundary of 
the region on which the space is dehned, i.e. s = 1 for arbitrary 9, oi 9 = n/2 for arbitrary s. 


21 













0.4 



Figure 9: The trajectory of figure 10 in (s, 6 ) plane. The horizontal axis is log(s) and the 
vertical axis is O/tt. The horizon corresponds to the left edge at log(s) = 0. 



Figure 10: Some generic solutions to the embedding equation of a D3’-brane in a finite 
temperature AdS^ x 5*5 background. The full embedding is cylindrically symmetric with 
respect to rotation around the Xg axis. The black disk represents the region behind the 
event horizon. 
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Figure 11: The trajectories of figure]^ in the {s,9) plane. The horizontal axis is log(s) and 
the vertical axis is O/tt. 


These solutions are referred to as the black hole and the Minkowski embeddings, respectively. 
These trajectories are specihed uniquely by giving the starting position of the trajectory 
along the boundary of the (s, 9) plane because the quadratic term in the equation of motion 
inferred from (2.15) degenerates there. 


The trajectory resulting from these initial conditions are illustrated in hgure 11 The 


same solutions in the original coordinates is illustrated in hgurej^ The trajectories 

starting at s = 1 boundary are the black hole embedding, and the trajectories starting at 
9 = 71/2 are the Minkowski embedings. 


B Retarded Green Function and Quasi-Normal Mode for 9{U) = 0 
background 


In this appendix, we will outline the computation of retarted Green function for the 9{U) 
fluctuation using the prescription of 


We start by generalizing (2.15) to include momen¬ 


tum and energy and write 


Sdbi = j7 A^ [ d^xdUUcos{9{U))J{l + fU^)9'{U)^+ ^^^^"^ . (B.l) 

{27rygs J y 

In order to compute the retarded Green function, we need to expand 9{U) to quadratic order 
around a solution to the equation of motion. This is quite complicated for generic solution 
9o{U), but is tractable for the trivial solution 9q{U) = 0. For that case, the linearized action 
reads 

Sun ~ \u [u^f9'{Uf - 9{Uf + (B.2) 

The equation of motion resulting from this action can be written in a canonical form in 
terms of 

9{x) = {x- lY^^y{x) , X = ^ (B.3) 
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with 


bJ = 


UsW 

yx’ 


k = 


Usq 


(B.4) 


so that the equation of interest is expressed as 

, / I 1 — iw/2 

0 = y'\x) < ' 


y\x) 


1 — X X 

k‘^ — w"^ + 2 — 2 q^ + w‘^ — 2 


8 (x - 1) 


8 x 


+ 




+ 


w 


8 (x + 1) 16a;2 


y(x) 


(B.5) 


For non-zero q^, this equation has four singular points at a; = 0, x = 1, a; = —1, and x = oo, 
and as such is not analytically tractable. But when q^ = 0, the singularity at x = — 1 goes 
away, and one arrives at a hypergeometric equation solved by 


, . . , _ tw _ tw I j- 

y{x) = ci(-l) ix 4 2 F 1 ( - 


fl 1 iw 

i+ir- 4 - 4 r+ 2 '^-T'" 


IW IW I J. 

-FC2(-1) 4 X 4 2 F 1 ( - 


1 I 
7 “ T I 

4 4 


1 i\ 1 iw 

4 + 4 r+ 2 'T + ^’" 


(B.6) 


One of the solution satisfies the infalling boundary condition at the horizon and the other is 
outgoing. We can then find the large U asymptotics the ingoing solution and find it to be 
of the form 

e{U) = A\og{U/Us) + B (B.7) 


where 
A = 

B = 


2^f2e^Y (1 - f) 


-V. 


(1 - f ) (V-l“> (i - ( J + I) d + ((i - i) + 5) + 7 + log(2)) 


(B.8) 




V. 


from which we infer that 


Gr{w) — — — log{Us/UcT) 


(B.9) 


with 


- = - f i - ( i + 1 I u. 


1 t 

A 2V" V2 V4'4;'~; " VV4 4 ;" 2/ / 


_^( 0 ) 


I \ 1 

T - T W + V 


+ 1] - 27 -log( 2 )^ . (B.IO) 


We see that the poles of the retarded Green function is encoded in the poles of 'ip^{x) 
function and is independent of Uqt, but the Green function itself is dependent on Uct 


through momentum independent contact terms as is shown in (B.9). 
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Holographic dictionary and the effective action for the order 
parameter 


In this article, the holographic dictionary (3.2) and (3.6) 


= Zbuik[C{x)] = j [D0((7, (C.l) 

played a critical role in providing an interpretation of brane embeddings in the bulk of space 
time in terms of field theory observables. The dependence on control parameter/boundary 
condition C{x) is somewhat implicit in the path integral expression on the right most side 


of (C.l). In this appendix, we will provide a formal path integral manipulation to make this 


explicit, as well as derive a formal path integral expression for the effective action for the 
field of the expectation value of the operator M{x) ~ {0{x)) corresponding to the bulk field 

e{U,x). 

The trick is to introduce an M{x) as an auxiliary field which when integrated out repor- 
duces the original path integral as follows. 


■'bulk 


[^(a;)] = [[D9{U,x)][DM{x)]e ^c{x)m{x))) 


where we are working in the path integral with cutoff at ?7 = Uuv fo regularize the contribu¬ 
tion of the counter-term although we will take the Uuv oo limit in the end. The auxiliary 
field lives only at f/ = Uuv- The presence of the holographic counter-terms guarantees that 
this limit is smooth. Because a logarithm is involved, we have intrduced yet another scale 
f/** although the dependence on this scale drops out in the Uuv oo limit. It is clear that 


integrating out M{x) will impose the boundary condition and reproduces (C.l). 


The expression (C.2) is useful for a variety of reasons. First, note that functionally 


differentiating with respect to C{x) pulls down a M{x). In that sense, we immediately 
associate M{x) with the expectation value {0{x)). 

We can also read off a path integral expression for the effective action of M[x) easily as 
follows. 

= J[D9{U, ^. (C.3) 

In this expression, the term linear in M{x) is imposing a boundary condition for the path 
integral over 9{U, x). It is also clear that T[M{x)] and W[C{x)] = — log(Z[C'(a;)]) are related 
by the standard Legendre transform at the leading order in saddle point approximation. 


whose correction can be computed systematically 35 
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The effective action r[M(a;)] is a complicated expression which includes terms with arbi¬ 
trary orders of M{x) and its derivatives. But it is formally defined umambiguously in (C.3) 
and can be computed systematically as an expansion in M{x) and its derivatives. 
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